Evidence for gap anisotropy in CaC6 from directional point-contact spectroscopy 
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We present the first results of directional point-contact spectroscopy in high quality CaC6 samples 
both along the ab plane and in the c-axis direction. The superconducting order parameter A(0), 
obtained by fitting the Andreev-reflection (AR) conductance curves at temperatures down to 400 
mK with the single-band 3D Blonder- Tinkham-Klapwijk model, presents two different distributions 
in the two directions of the main current injection, peaked at f .35 and 1.71 meV, respectively. By 
ab-initio calculations of the AR conductance spectra, we show that the experimental results are in 
good agreement with the recent predictions of gap anisotropy in Ca,C%. 

PACS numbers: 74.50.+r, 74.45.+C, 74.70.Ad 



The discovery of a relatively "high- IV' superconduc- 
tivity in graphite intercalated with Ca [HQ, Yb [l| and, 
very recently, Sr [3j, |4[ has strongly revived the inter- 
est in the Graphite Intercalated Compounds (GICs) and 
their physics. The Ca-intercalated graphite, CaCg, with 
its "record" T c of about 11.5 K, in particular, has been 
the subject of various theoretical and experimental in- 
vestigations in the past two years (for a short review of 
the initial results see Q). One of the most important 
questions, however, is still not clear: what is the mag- 
nitude and anisotropy of its superconducting gap? The 
first experiments (STM, penetration depth, specific heat) 
on CaCg have evidenced a single, apparently isotropic, s- 
wave gap with a ratio 2A/kftT c of the order of the BCS 
value laM 0- Recent tunnel spectroscopy results, on the 
other hand, claimed the presence of an isotropic gap with 
a magnitude more than 40% higher than that reported 
earlier [§] . The spread of gap values measured up to now 
range between 1.6 meV Jg| and 2.3 meV The im- 
portant point is that all these experiments have either 
probed a bulk property |a| or a directional one alon 
the c-axis direction @, 0,111 ■ As pointed out in Ref 
the presence of anisotropic or two-gap superconductivity 
in CaCe cannot be ruled out until tunneling or point- 
contact measurements are performed also along the ab 
direction. On the other hand, recent first-principles den- 
sity functional calculations of the superconducting prop- 
erties of CaC6 have supported the presence of a moder- 
ately anisotropic gap which varies between 1.1 and 2.3 
meV, depending on the k-point and the ir or interlayer 
(IL) sheet of the Fermi surface (FS) involved [l(|. Such 
an anisotropy can be revealed by directional spectroscopy 
measurements performed along both c and ab direction. 

In this paper we present the results of point-contact 
Andreev-reflection (PCAR) spectroscopy performed on 
high-quality bulk samples of CaC6 0. By using a spe- 



cial technique to realize the contacts, that proved very 
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successful and effective in the case of MgB 2 
were able to perform directional PCAR spectroscopy at 
very low temperature both along the ab plane and the 
c-axis direction. Two different gap distributions in the 
two directions can reproducibly be extracted from the ex- 
perimental data. When compared to the results of new 
first-principles calculations these findings unequivocally 
prove the anisotropy of the superconducting gap in CaCg. 

The high-quality CaCg bulk samples used for our mea- 
surements were synthesized by reacting highly oriented 
pyrolytic graphite (with a spread of the c axis orienta- 
tion < 0.4°) for several weeks at 350°C with a molten 
alloy of Li and Ca Q . The resulting CaCg samples have 
a shiny golden surface. They are very sensitive to air and 
moisture which rapidly damage the sample surfaces. X- 
ray analysis has shown mainly the CaC6 reflections with 
a small (< 5%) contribution from impurity phases. Fur- 
ther details on the characterization of the samples may be 
found in Ref. [||. All samples used for PCAR spectroscopy 
(size w 1 x 1 x 0.2 mm 3 ) were selected to have a very sharp 
superconducting transition (AT C (10% - 90%) = 0.1 K) 
with the onset at T c = 11.4 K. 

The point contacts were made by using a non- 
conventional technique we called "soft" PCAR spec- 
troscopy 11 , 13] ■ Instead of using the standard metallic 
tip, a very small (0 ~ 50 /im) drop of Ag conductive 
paint, put on the etched or freshly cleaved surfaces of 
the sample is used as a counterelectrode. Such contacts 
are particularly stable both in time and towards tem- 
perature variations and they allow to inject the current 
mainly perpendicular to the contact plane. A fine tuning 
of the junction characteristics at low temperature can be 
done by applying short voltage or current pulses. Fur- 



ther details on the technique can be found in Refs. 1 1 lL 1 1 3 . 
Due to the mentioned high sensitivity of CaCg samples' 
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FIG. 1: (Color online) (a) Raw point-contact conductance 
curves of various aft-plane contacts at 4.2 K. For clarity the 
curves are vertically shifted of the amount shown close to each 
curve, (b) The same as in (a) but for various contacts with 
current injection mainly along the c axis. In each panel a 
sketch of the contact geometry is also shown. 

surface to air, the room-temperature preparation of the 
contact was done inside a sealed glove bag filled with 
pure He gas or in a glove box with Ar atmosphere. Af- 
ter the contact was made the junction was very rapidly 
transferred to the cryostat in a sealed container. Con- 
tacts were made either on the flat aft-plane surface or on 
the narrow lateral side of the samples. Referring to the 
main direction of current injection, we call them c-axis 
and aft-plane contacts, respectively (see insets of Fig. 1). 

The conductance curves, dl/dV vs. V, were obtained 
by numerical differentiation of the measured I — V curves 
and subsequently normalized by dividing them by the 
normal-state conductance measured at T > T c . For this 
reason, in all the contacts, we therefore carefully studied 
the temperature dependence of the conductance in order 
to determine the critical temperature of the junction, i.e. 
the 'Andreev critical temperature', T C A . In an overall of 
35 different contacts, T A was found to be 11.3 ±0.1 K, in 
best agreement with the bulk T c 's of the samples and in 
contrast with a previous report 0]. This fact proves the 
high quality of samples and surfaces in the contact re- 
gion. For simplicity, we will therefore refer to the critical 
temperatures of the contacts as T c in the following. 

Fig. Q] shows several raw conductance curves as func- 
tion of bias voltage measured both in aft-plane contacts 
(a) and in c-axis ones (b) at 4.2 K. The curves show clear 
Andreev-reflection (AR) features, an almost flat conduc- 
tance (at V > 8 — 10 meV) and no dips that usually 
are a sign of the failure in reaching the conditions for 
pure ballistic conduction in the contact [3, EH]. The 
normal resistance i?N of all the good contacts is between 
0.75 and 6.4 fi. By knowing the mean free paths and 
the residual resistivities of CaCg along the aft plane and 
in the c-axis direction, i.e. £ a b = 74 nm, l c = 4.7 nm, 



Po.ab = 0.8 jufi-cm and po,c = 24 /iSl-cm [16J, [17J we can 
apply the Sharvin formula for the contact resistance in 
the ballistic regime in order to determine the contact 
radius a = (4p ^/37ri? N ) - 5 0. The condition for full 
ballistic transport (a <C £) is totally verified in aft-plane 



contacts, where a a b ~ 6 — 18 nm. In c-axis junctions, 
where a c « 14 — 24 nm but the conductance curves do 
not show any sign of heating [l3| , the presence of at least 
30 parallel contacts in the junction area is expected. 

After normalization, the conductance curves were 
fitted to the modified 3D Blonder-Tinkham-Klapwijk 
(BTK) model 0, 0, H(J . In the single-band form it con- 
tains three fitting parameters: The gap A, the barrier- 
height parameter Z and the broadening T which accounts 
for both intrinsic (quasiparticle lifetime) and extrinsic 
phenomena that broaden the AR conductance [l|| • 

In order to increase the experimental resolution of our 
measurements we decided to perform part of the PCAR 
experiments at very low temperature in a Quantum De- 
sign measurement system (PPMS) with 3 He insert. 

Fig. [2] (a) shows the normalized conductance curves 
(circles) of a typical aft-plane contact at various temper- 
atures from 400 mK up to T c . At any temperature the 
single-band 3D BTK model fits the data very well (solid 
lines). At the lowest T, the values of the fitting param- 
eters are: A = 1.44 meV, T = 0.61 meV and Z = 0.75. 
In panel (c) we display the order parameter A obtained 
from the data given in (a). Its temperature dependence 
almost perfectly follows the BCS-like expression (solid 
line) with 2A(0)//cbX , c = 2.98 which is sensibly smaller 
than expected from BCS theory. 

In Fig. [2] (b) and (d) we report the same data for a 
c-axis contact. As for the aft-plane case, the curves are 
well fitted by the single-band 3D BTK model which gives 
at 400 mK: A = 1.7 meV, T = 0.84 meV and Z = 0.97. 
The temperature dependence of A is very close to the 
expected BCS one with a ratio 2A(0)/fc B T'c = 3.48, in 
best agreement with the weak-coupling BCS value. 

It is worth noticing that the Z values observed in c- 
axis contacts (between 0.74 and 1.01) are systematically 
greater than those of aft-plane contacts (between 0.48 
and 0.75). According to the 3D BTK model [H, this 
difference can be explained by the different Fermi veloc- 
ities of CaC6 in the aft plane (w a b = 0.54 x 10 6 m/s) and 
along c axis (v c = 0.29 x 10 6 m/s), thus confirming the 
directionality of our point contacts. 

The AR curves shown in Fig. [1] and [2] are rather 
small in amplitude, as already observed in all the "soft" 
PCAR measurements on MgB2 and related compounds 
11, [H, resulting in T values substantially greater 
than those expected for the quasiparticle lifetime. As re- 
cently observed in lithographically fabricated Cu-Pt-Pb 
this additional broadening can be ex- 
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nanocontacts 

plained by the presence of pair-breaking effects induced 
by the scattering in a thin disordered layer present at the 
NS interface. This is the case of our point contacts, due 
to a disordered layer on the surface of Ag grains that also 
makes the residual resistivity of the paint be five orders 
of magnitude greater than in pure Ag. 

The reproducibility of the PCAR data was very good. 
Most of the contacts, obtained both in 4 He and in 3 He 
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FIG. 2: (Color online) Normalized dl/dV vs. V curves at different temperatures down to 400 mK in an afc-plane contact (a) 
and in a c-axis one (b) (open circles). Solid lines: best-fit curves given by the single-band 3D BTK model. Panels (c) and 
(d) show the temperature dependency of the order parameter A (full circles) in the afe-plane direction and in the c-axis one, 
respectively, as determined from the BTK fits shown in (a) and (b). Solid lines are the BCS-like fits. 



cryostat, show dl/dV curves and temperature dependen- 
cies of quality similar to that presented in Fig. [2] In 
15 ao-plane contacts the order parameter A(0) ranged 
between 1.1 meV and 1.7 meV with the distribution 
shown in Fig. [3] (a). In 14 c-axis contacts A(0) ranged 
between 1.3 meV and 1.94 meV with the distribution 
shown in Fig. [3] (b). The figure also shows the Gaussian 
curves that best fit the distributions. They are peaked at 
A a b(0)=1.35 meV and A c (0)=1.71 meV and show stan- 
dard deviations s a b=0.14 meV and s c =0.08 meV, respec- 
tively. The results in c-axis direction are in very good 
agreement with the gap values previously reported in 
Ref. @, 0- A minority of contacts (3 in aft-plane and 
3 in c-axis direction) have shown gap values between 2.1 
and 2.4 meV, similarly to the results of Ref. [tj 

The complex microscopical nature of our point con- 
tacts leaves some uncertainty about the true direction 
of current injection, particularly in the case of contacts 
on the side faces of the sample (i.e. a&-plane contacts) 
where, due to the intrinsic inhomogeneity of the cleaved 
surface, current injection along c axis is also possible. 
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FIG. 3: (Color online) Distributions of the different A(0) val- 
ues measured in afc-plane contacts (a) and in c-axis ones (b) 
at 4.2 K (red) and at 400 mK (light red). Dash black lines 
are the fits of the total distribution to a Gaussian curve. 



However, the clear difference observed between the most 
probable A(0) values in a6-plane and c-axis contacts pro- 
vides strong evidence for a gap anisotropy in CaCg. 

In order to compare our results with the theoretical 
predictions of gap anisotropy in CaC6 [ic} we calcu- 
lated the Andreev-reflection conductance curves by first- 
principles methods. We have a SN junction where S = 
CaC6 and N = Ag. Let's label with the suffix i — 1,2,3 
the three sheets of the CaCg Fermi surface (FS) (ir and 
interlayer (IL) bands). If n is the unitary vector in the 
direction of the injected current, i>ik, n = Vjk • n are the 
corresponding components of the Fermi velocities in the 
superconductor at wave vector k for band i-th. Tak- 
ing into account that Ag has a quasi-spherical FS and 
an almost constant Fermi velocity i>n ^ «ik, the corre- 
sponding quantity in the normal metal will be i>N,n = i>n- 
Following Rcfs.l22l.l23l we finally obtain the total AR con- 
ductance as: 



Ei(^ikn(£)- 



lt[»ik,o+«N] 



i'Si 



(1) 



/FSi 



where: Ofsi is the integral over the i-th FS, i.e. 

( , Vik - n l nr )FSi = f ^ n f % V2 S(E ik )d 3 k. 

Cikn(-E) is the BTK conductance of the i-th band ex- 
pressed in terms of the quantities N^(E) = E/^JE 2 — A? k 
and Nfr(E) = A ik /\/i5 2 — A? k (whose real parts are the 
quasiparticle and the pair density of states in the same 
band, respectively) and of Z n values. Aik is the gap value 
for band i-th at point k over the FS, recently calculated 
from first principles [Io| . The values of Z n used in the 
calculation are taken similar to those of the curves shown 
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FIG. 4: (Color online) Theoretical AR conductances calcu- 
lated at T — by Eq. (1). (a) Current injected along the 
a axis with Z = 0.75 and T = (black) and V = 0.6 (red); 
(b) current injected along the c axis with Z = 1 and V = 
(black) and F = 0.8 (red). Experimental curves at 400 mK 
are shown for comparison (blue circles). 

in Fig. [21 i.e. Z a b = 0.75 and Z c = 1. The explicit ex- 
pression of (Jikn(-B) can be found in Ref. |2fJ. 

The results of the calculations given by Eq. (1) are 
shown at T = K and for a- and c-axis directions (top 
curves in Fig. [4] (a) and (b)). The conductance calcu- 
lated along the b direction is almost identical to the one 
in a direction. At T = the topology of the CaC6 FS 
and the calculated anisotropy of the 7r and IL gaps re- 
sult in a sizeable anisotropy of the AR conductance. In 
ab direction, it exhibits a sharp peak (related to the it 
gap) at about 1.38 meV and a broad shoulder (mainly 
related to the IL gap) at about 1.9 meV. In c direction, 
as expected from the shape of the FS, the role of the IL 
(Ca) gap becomes more important and the conductance 
shows two distinct peaks with almost the same height. 
However, for T ^ even at T = K these anisotropic 
features are rapidly smeared out. The middle curves of 
Fig. 0Ja) and (b) show the effect on the theoretical con- 
ductances of a broadening similar to that observed at 
very low T in the experimental curves of Fig.O The con- 
ductances become similar to single-gap ones and can be 
perfectly fitted by single-gap 3D BTK curves, as shown 
in Fig. |4] (black dash lines). The use of more complex 
fitting models (anisotropic or two-band BTK) that we 
tested on our data does not improve the fit substantially, 
as already pointed out in Ref. S 0. In Fig. 1 the ex- 
perimental a&-plane and c-axis conductances measured 
at 400 mK are included too (circles) in order to show 
the remarkable agreement with the theoretical curves for 
the same level of broadening. Although this broadening 
washes out the fine anisotropic structures of the conduc- 
tance, a clear sign of the underlying gap anisotropy is 
still present since the 3D BTK fit of the theoretical con- 
ductances gives different order parameters in the two di- 
rections, A = 1.5 meV (r = 0.65 meV, Z = 0.765) for 
a-axis direction and A = 1.7 meV (r = 0.92 meV, Z = 
1.015) for c-axis one. The c-axis value is in perfect agree- 
ment with the experimental results (both single curves at 



400 mK and the peak of the distribution of the 14 dif- 
ferent contacts). In the a6-plane case, the experimental 
A values from the curves at 400 mK and from the peak 
of the distribution of Fig. [3] (a) (ranging from 1.3 meV 
to 1.44 meV) are smaller than the value obtained from 
the fit of the theoretical conductance. This discrepancy 
could be ascribed to a possible slight overestimation of 
the small n gap (and, maybe, an underestimation of the 
large IL gap associated with Ca FS) in the theoretical 
calculations. This fact appears reasonable if one consid- 
ers that the first-principle calculations of Ref. [lj| led to 
an underestimation of T c of about 17 %. 

In conclusion, the first directional PCAR measure- 
ments in CaCg carried out also at T = 400 mK both 
along the a&-plane and the c-axis direction give strong 
and reproducible evidence of the predicted anisotropic 
nature of the superconducting gap in this GIC. New first- 
principles calculations of the expected anisotropy in the 
AR conductance curves fully support this conclusion and 
indicate that the actual gap anisotropy in CaC^ could be 
even slightly greater than theoretically predicted. 

We thank Lilia Boeri, O.V. Dolgov, E.K.U. Gross and 
LI. Mazin for useful discussions. This work was done 
within the projects: PRIN (No. 2006021741) and Cy- 
bersar (cofunded by MUR under PON). V.A.S. acknowl- 
edges the support of RFBR (Proj. No. 06-02-16490). 



[1] T. E. Weller et ai, Nature Phys. 1, 39 (2005). 
[2] N. Emery et ai, Phys. Rev. Lett. 95, 087003 (2005). 
[3] J. S. Kim et ai, Phys. Rev. Lett. 99, 027001 (2007). 
[4] M. Calandra and F. Mauri, Phys. Rev. B 74, 094507 
(2006). 

[5] I. I. Mazin et ai, Physica C 460-462, 116 (2007). 

[6] N. Bergeal et ai, Phys. Rev. Lett. 97, 077003 (2006). 

[7] G. Lamura et ai, Phys. Rev. Lett. 96, 107008 (2006). 

[8] J. S. Kim et ai, Phys. Rev. Lett. 96, 217002 (2006). 

[9] C. Kurter et ai, Phys. Rev. B 76, 220502 (2007). 
[10] A. Sanna et ai, Phys. Rev. B 75, 020511(R) (2007). 
[11] R.S. Gonnelli et ai, Phys. Rev. Lett. 89, 247004 (2002). 
[12] R.S. Gonnelli et ai, Phys. Rev. Lett. 97, 037001 (2006). 
[13] D. Daghero et ai, Phys. Rev. B 74, 174519 (2006). 
[14] Goutam Sheet, S. Mukhopadhyay, and P. Raychaudhuri, 

Phys. Rev. B 69, 134507 (2004). 
[15] Y.G. Naidyuk, I.K. Yanson, Point-Contact Spectroscopy, 

Springer Series in Solid-State Sciences, Vol. 145, 2004. 
[16] J. S. Kim, unpublished results. 

[17] A. Gauzzi et ai, Phys. Rev. Lett. 98, 067002 (2007). 
[18] G.E. Blonder, M. Tinkham and T.M. Klapwijk, Phys. 

Rev. B 25, 4515 (1982). 
[19] A. Plecenik et ai, Phys. Rev. B 49, 10016 (1994). 
[20] S. Kashiwaya et ai, Phys. Rev. B 53, 2667 (1996). 
[21] P. Chalsani et ai, Phys. Rev. B 75, 094417 (2007). 
[22] I.I Mazin, Phys. Rev. Lett. 83, 1427 (1999). 
[23] A. Brinkman et ai, Phys. Rev. B 65, 180517(R) (2002). 



